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Abstract: Since Global Navigation Satellite Systems (GNSS) show degraded perfor- 
mance in dense urban and indoor areas, a positioning sensor based on Digital Video 
Broadcast-Terrestrial (DVB-T) systems is presented in this paper. DVB-T signals can be 
considered as signals-of-opportunity for positioning, due to their good properties. One of the 
challenges to overcome is to distinguish the signals from different emitters. Here, we suppose 
that the user can first compute his position by GNSS during an initialization phase, which is 
used for solving all the ambiguities concerning DVB-T emitters. Starting from there, DVB-T 
signals can be used for aiding positioning when the user enters a GNSS-blocked area, up to 
a limit case, where all the GNSS satellites are not in view and only DVB-T signals are used 
for positioning. We tested this method by simulation, by adopting the Hata model for the 
emitter attenuations and the Rayleigh model for multipath. The obtained results show good 
performance if the receiver correctly associates the signal to the user's motion. 

Keywords: DVB-T; signal-of-opportunity; GNSS; positioning 



1. Introduction 

Global Navigation Satellite Systems (GNSSs) have been widely used in many applications for 
positioning, navigation and timing. However, in urban areas and indoor environments, where position 
information is required for many applications, GNSSs show degraded performance in terms of precision 
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and availability, because of the signal loss or attenuation and multipaths due to obstacles. Fortunately, 
in these areas, many local networks are deployed, such as 2G, 3G, WiFi, LTE (Long Term Evolution) 
and DVB-T (Digital Video Broadcast-Terrestrial). They were originally designed for other purposes, but 
they can be used for positioning, thanks to their properties, such as a high signal-to-noise ratio (SNR). 

In this paper, we consider DVB-T systems. DVB-T is the European digital TV standard. It adopts 
the Orthogonal Frequency Division Multiplexing (OFDM) technique. Therefore the whole bandwidth is 
divided into many subcarriers, in which the pilot subcarriers are included. We used these pilot subcarriers 
to estimate the ranges between the receiver and different emitters with a mechanism similar to the one 
used by GNSS receivers. Since the SNR requirement for ranging is much lower than the one required 
by TV service, the receiver is able to see several emitters in one point. If three or more signals are 
successfully processed, the receiver can provide a DVB-T only the positioning; otherwise, it can be used 
to assist GNSS. 

According to the document [1], two different network types are introduced, Multi-Frequency 
Network (MFN) and Single Frequency Network (SFN). In MFN, different emitters transmit the same 
signal on different frequencies, and they are not exactly synchronized; while in SFN, all the emitters 
simultaneously transmit the same signal in the same frequency. The synchronization is provided by 
some professional GPS timing receivers, and this represents a key point when DVB-T signals are used 
for positioning purposes. Notice that the transmission of the emitter identifier (ID) is optional: this is 
not an issue for the provision of the DVB-T service, but it is a problem to be solved when the signals are 
used for positioning. 

It is interesting to note that some companies transmit the emitter ID within their bit streams, which 
are then completely equal apart from this small difference (this is done, for example, by Rai in Italy). 
This causes a (very limited) penalty, but is very useful for network management and control. It is clear 
that this extra-information can be very useful also for positioning purposes, since it highly simplifies 
the association of each echo to the corresponding emitter. Anyway, since the ID emitter transmission is 
optional in the DVB-T standard (which is much more used and available than DVB-Handheld (DVB-H)), 
we have not considered it in the study, and we have proposed another way to solve the ambiguity. 
(Clearly, some actions could be adopted in the future to convince the TV companies to transmit it, to 
further exploit DVB-T for positioning.) 

The scenario analyzed in this paper considers a position device (PD) consisting of a hybrid 
GNSS/DVB receiver, where GNSS is the primary positioning system, and DVB-T is used as a back-up 
when the number of GNSS satellites is not sufficient for computing the receiver position. In the first 
phase (called the initialization phase in the rest of the paper), we suppose that the PD is able to compute 
its position by using GNSS only. In this phase, the PD keeps sensing the DVB-T spectrum in order 
to identify all the DVB-T emitters and computes and tracks the ranges between each emitter and the 
PD itself. Since the emitter positions are known, the ranges can be used to correctly associate each 
received signal to the corresponding emitter. Notice that a problem could arise when two or more ranges 
assume the same value; however, the probability that such an event continuously persists as long as the 
PD is in open sky is very low and can be neglected. When the PD enters a GNSS-hostile zone (indoor 
environment, urban or natural canyon or similar) and the number of visible satellites becomes lower 
than four, the PD can use the DVB-T ranges to integrate them, up to a limit solution, where no satellites 
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are visible and the position is computed by using DVB-T data only. This last case obviously represents 
an extreme situation, which may become very fragile, in the presence of path ambiguities encountered 
when two ranges become equal. The problem can be quite completely eliminated if some redundancy 
is available or when auxiliary systems are available (i.e., map matching). In Section 4, two methods are 
presented to solve the path ambiguity problem, with the assumption that one GNSS satellite is visible. 

In the literature, the problem of using DVB signals for pseudo-range calculation is addressed in some 
papers. In [2], a pseudo-range measurement method in a single channel based on the DVB-T signal is 
presented. The method uses the classical functions of an OFDM receiver to demodulate the OFDM signal 
and estimate the Carrier Frequency Offset (CFO) and the Sampling Clock Offset (SCO). The method has 
been assessed in [3,4] through experimental results using the French digital TV terrestrial network based 
on the DVB-T standard. In [5], a pseudo-range measurement method based on demodulated data is 
described. However, no positioning algorithm is presented, because of the fact that all the emitters in an 
SFN simultaneously transmit the same signal on the same frequency, and the problem of distinguishing 
the signals from different emitters is not addressed in the paper. A method to distinguish the emitters 
is proposed in [6,7], based on the idea of introducing some artificial delays on the system side of a 
DVB-SH (satellite services to handhelds) system. Since DVB-H and DVB-T have a similar structure, 
in principle, the method could be also studied for DVB-T, but since this option is not available in the 
current DVB-T standard, the method has been ignored in our study. In our previous works [8-10], a 
method to distinguish different signals in a dynamic scenario has been studied by assuming that the 
PD is a hybrid DVB/GNSS receiver working in DVB-mode after an initialization phase provided by 
GNSS. In this phase, the initial position and velocity are available and represent the key point to allow 
the trilateration operations based on DVB-T signals. In [8], DVB-T is used as an assistant of GNSS; no 
DVB-T-only solution is provided. In [9], only an ideal scenario with no signal attenuation and multipath 
has been considered. In [10], only a simple channel model is used, and no path ambiguity solutions 
are considered. In this paper, these aspects are taken into account, by adopting the Hata model [11] 
for the attenuation and a Rayleigh model for the multipath. Two different path ambiguity solutions are 
also explained in this paper. Moreover, the impact of unsynchronization of different emitters is also 
investigated. The performance in terms of position errors in different SNR is analyzed. Some other 
minor, but important, points are also investigated in this paper, such as Doppler effects and a windowing 
technique to mitigate the side lobes. 

This paper is organized as follows. Section 2 presents a description of SFN DVB-T. Section 3 
describes the pseudo-range estimation method. Section 4 presents two methods to solve the path 
ambiguity problem. Section 5 presents the simulation results. The paper ends with the conclusions 
and a future work description. 

2. SFN DVB-T Description 

According to the DVB-T document [1], DVB-T systems adopt the OFDM technique. Four different 
types of subcarriers are used, among which the pilot subcarriers exist. They are used for channel 
estimation in DVB-T applications, but they can be also used for positioning by correlating the incoming 
signal with a local-generated replica, with a method similar to that used in GNSS receivers. In order to 
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avoid inter-symbol interference, the Cyclic Prex (CP) is introduced in the OFDM system. This property 
is very helpful for our positioning purpose. In the coming part, we will give a brief introduction of 
the SFN DVB-T system, highlighting the points that impact more on the capabilities of these signals to 
behave as Signals of Opportunities (SoOs) for positioning. 

2.1. OFDM 

OFDM is a digital multi-carrier modulation method. It divides the bandwidth into a large 
number of closely-spaced sub-bands. On each sub-band, a Quadrature Amplitude Modulation 
(QAM) or Quadrature Phase-Shift Keying (QPSK) is used, and all the sub-band signals are 
summed together. The symbol rate is the same for each subcarrier and is equal to the sub-band 
bandwidth. This way, the subcarriers are orthogonal, and the spectra can overlap without causing 
Inter-Carrier Interference (ICI) when the receiver is well synchronized. To avoid Inter-Symbol 
Interference (ISI) in multipath fading channels, a guard interval is inserted prior to the OFDM symbol. 
This interval is used to transmit an exact replica at the end of the OFDM symbol, referred to as CP. This 
gives the OFDM system an excellent multipath resistance: the receiver can easily avoid the ISI if the 
multipath time- spreading is shorter than the guard interval. 

The Fast Fourier Transform (FFT) algorithm can be used to implement the OFDM modulation, and 
the good efficiency of this algorithm allows a large number of subcarriers in operation. In Figure 1, the 
block diagram of an OFDM transmission system is presented. The output of the parallel to the serial 
block can be expressed as: 

NppT—l 

= IFFT(c£) [n] = — 4^>[.i' 2 *- 

with 0 < n < Nfft — 1 

where: 

• k is the OFDM symbol index, 

• n is the sample index, 

• p is the subcarrier index, 

• N FFT is the FFT size, 

• sjj is the n-th sample of the £>th OFDM symbol and 

• Cp is the data (constellation complex value) on the p-th subcarrier of the A;-th OFDM symbol. 

Normally, in OFDM modulation, the number of subcarriers, N FFT , is chosen as a power of two in 
order to efficiently use the FFT algorithm. Before the samples are transmitted, the CP is inserted at the 
beginning of each OFDM symbol. Therefore, the emitted signal with the CP is: 

„k = ~k 

°n °rnod{n+NpFT~ Nqp ,Npft) ^2"^ 

with 0 < n < N FFT + N C p - 1 

where: 

• is the n-th emitted sample of the k-ih. OFDM symbol and 



1\ FFT \ i\fft/ (1) 
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N C p is the length of the cyclic prefix guard interval. 



Figure 1. Orthogonal Frequency Division Multiplexing (OFDM) transmission system. 
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On the receiver side, the CP has to be removed. To do this, a coarsely-timed synchronization is needed 
to find the starting instant of the OFDM symbol. The synchronization can be achieved, for example, by 
using the Van de Beek algorithm [12]. This algorithm is based on the fact that the guard interval is 
a replica of the end of the OFDM symbol. The correlation of two groups of N cp samples spaced by 
N FFT is calculated, and the absolute maximum of the correlation is used to estimate the beginning of 
the OFDM symbol. The synchronization can be achieved more precisely by adopting a delay estimation 
method. If the starting point of the FFT window is slightly wrong, but still within the CP, the receiver still 
has all wanted samples. These samples are circularly shifted, comparing the samples originally emitted. 
This cyclic permutation only affects the phase of the received symbol output by the FFT. This additional 
phase rotation is typically compensated for by the frequency equalizer, since it is summed to the phase 
rotation introduced by the channel. 



2.2. DVB-T System 

DVB-T is a digital broadcasting standard created by the European Telecommunications Standards 
Institute. For mobile services, another available standard is DVB-H (Handheld), derived from DVB-T. 
In this paper, only the DVB-T is considered, but the proposed positioning technique can be also used for 
the DVB-H, which is more suitable for mobile users in dynamic scenarios. 

The DVB-T standard family has adopted the OFDM modulation to provide high data rates along 
with robustness against multipath. The transmitted signal contains four types of subcarriers, as 
described hereafter: 

• Null subcarriers, which are placed on the edge of the signal spectra and are not used. They serve 
as frequency guard bands to avoid the out-of-band emission of the OFDM signal. 

• Data subcarriers, which are the payload of the OFDM symbol. Each subcarrier carries some data 
bits, depending on the modulation scheme. The data subcarriers in one OFDM frame of a DVB-T 
are modulated using either QPSK or 16-QAM or 64-QAM. 
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• Transmission Parameter Signaling (TPS) subcarriers, which are used for transmitting signaling 
parameters related to the transmission scheme. These subcarriers are Binary Phase-Shift Keying 
(BPSK) -modulated at the normal power level. The TPS subcarriers may contain the cell-ID 
information, which is the only difference between signals transmitted from different emitters 
in SFN. This information can be used to associate the signals to the corresponding emitters. 
Unfortunately, this cell-ID information is optional; the transmitters can discard this identifier by 
setting them to zero. 

• Pilot subcarriers, which are used for channel estimation and equalization. They can be divided into 
two groups: the continuous pilots and the scattered pilots. The continuous pilots are always placed 
on the same subcarriers in all OFDM symbols. On the contrary, the scattered pilots are placed 
on different subcarriers in successive OFDM symbols. They are inserted every 12 subcarriers, 
and the first pilot is placed on one of four different subcarriers location (three, six, nine or 12), 
depending on the index of the OFDM symbol. All these pilot subcarriers are modulated by a known 
Pseudo-Random Binary Sequence (PRBS) with a 4/3 boosted signal amplitude compared to 
data and TPS subcarriers. The organization of the pilot subcarriers is shown in Figure 2. The 
positioning system described in this paper is based on the processing of the scattered pilot 
subcarriers (SPS). 



Figure 2. Digital Video Broadcast-Terrestrial (DVB-T) pilot organization. 
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In the process of positioning computation, we are interested in the ranges between the receiver and 
different emitters. These can be obtained by estimating the propagation delay. In DVB-T, the propagation 
delay is mainly related to the FFT size, N FFT (which corresponds to the number of subcarriers of 
the OFDM symbol), the ratio between the Cyclic Prefix length and the useful OFDM symbol length, 
Pcp = N C p/N FFT , and the sampling period, T samp . They are summarized in Table 1 according to the 
document, [1]. 

These values determine the maximum delay that a DVB-T receiver can handle. For example, in 
the 2 K mode with p C p = 1/4 and an 8 MHz channel, the maximum detectable delay is 56 [is, 
which corresponds to a path difference of about 16.8 km, assuming the propagation velocity is equal 
to the speed of light in vacuum (c = 299,792,458 m/s). Therefore, in general, it is possible to 
associate to each possible system configuration of Table 1 a maximum allowable path difference, Ap max , 
to avoid ISI. 
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Table 1. DVB-T parameters. 



Parameter Possible Values 



NppT 


2,048 (Mode 2 K) 




4,096 (Mode 4 K) 




8,192(Mode 8 K) 


PCP 


1/32, 1/16, 1/8, 1/4 


T S amp 


7/64 us (8 MHz)) 




1/8 us (7 MHz) 




7/48 us (6 MHz) 




7/40 us (5 MHz) 



As shown in Table 1, in the 2 K model with an 8 MHz bandwidth, the sampling interval is 
7/64 u.s. Multiplying by the speed of light, we have 32 m. By half, it is 16 m. This accuracy can 
be improved by some interpolation methods. Furthermore, we may have four or more DVB-T emitters 
in a 2D positioning; the obtained position resolution is acceptable for the urban scenario. 

3. Positioning Based on DVB-T 

It is well known that the position of an object capable of measuring the distances between itself and 
some reference points can be performed by using the trilateration method; this is the technique typically 
employed in GNSS [13], and the same method can be utilized when the references points are the emitters 
of DVB-T signals. To implement the method, the PD has to: 

1 . identify a reference frame, 

2. identify the emitter locations, pj = (xj, y-i, Zj), where i = 1, • • • , N e and N e are the number of 
visible DVB-T emitters, 

3. measure the ranges, r i; between each emitter and the PD 

4. write and solve the navigation equations. 

The implementation of the first and fourth tasks does not present any difference with respect to a 
classical GNSS receiver; so, it is not described in this paper. The second task is the most challenging, 
since the signals transmitted by the DVB-T emitters do not contain the station identifier; so, the receiver 
is not able to associate each received signal to a specific station. This is not an issue for TV reception, 
but it is a problem for positioning. Notice that the emitter identifier is foreseen by the DVB-T standard, 
but its transmission is optional; so, it is ignored in this study. The method proposed in this paper to solve 
this problem is described in Section 3.1. The third task requires the processing of an OFDM signal to 
estimate the range; so this has to be re-designed with respect to a GNSS receiver. 

The next sections will be devoted to the methods used to implement tasks 2 and 3, while the position 
computation can be implemented by solving the navigation equations with a classical extended Kalman 
filter (EKF) technique [14,15]. The experimental results shown in this paper have been obtained by using 
this approach. 
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3.1. Initialization Phase for DVB Ranging 

We assume here that in the future, a PD is likely to be equipped with a GNSS sensor as a primary tool 
and with other sensors able to exploit the nearby signals-of-opportunity. This means that the architecture 
of a PD will include a hybrid receiver, which computes its position with GNSS and resorts to SoSs only 
when the GNSS satellites are not visible. Moreover, in the case of DVB-T towers, it is realistic that the 
PD has a map of their locations (they are indeed available and easy to find also on the web). 

In this scenario, we can suppose that the trilateration with DVB-T-based measurements is generally 
preceded by a phase of position computation based on GNSS, performed when the PD is in open sky 
with complete visibility of GNSS satellites. This phase is the standard mode of operation and, at the 
same time, represents the initialization phase for the position computation based on DVB-T. During this 
phase, the PD continuously evaluates the ranges, r i5 between each emitter and the PD, by using its own 
position, provided by GNSS, and the location of the nearby DVB-T emitters, whose map is stored in the 
PD memory. Therefore, in this phase, the PD can also associate the estimated ranges with the emitter 
locations, solving the problem of emitter identification, as shown in Figure 3. 

Figure 3. Initialization pahse. 
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Each range can be written at each discrete-time instant, n, as: 

r i)G [n] = N i)G [n]r s + 5 i)G [n] (3) 

where N i:G [n] is an integer, r s is the distance traveled by the DVB-T signal during an OFDM symbol 
time, £i,gM is a fractional range with respect to r s and the subscript, G, stands for GNSS. Notice that 
Ni )G [n] could be also zero, depending on the coverage area of the DVB-T towers. At the same time, 
the PD can estimate the ranges from the DVB-T emitters by using the SoOs. This can be obtained 
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by measuring the fractional range, (^[n], as shown in Section 3.2, and by evaluating N^ G [n} from 
Equation (3), as: 



N i>G [n) 



(4) 



where [_-J stands for the integer part. Therefore, the PD knows at each instant, n, both r i G [n] and the 
range provided by the DVB-T SoOs, which can be written as: 



r itD [n] = N itG [n]r s + 5 i)D [n] 



(5) 



and can keep these values aligned. At this point, if the PD enters a GNSS-blocked area at the time 
instant, n 0 , it can start the position computation in DVB-mode by using the range measured at the time 
instant, n 0 — 1, that is, the range measured during the initialization phase. This represents a hot start for 
the PD, similar to the one encountered in a PD only working in GNSS-mode. 

Another benefit of the initialization phase is that the out bound peaks can be discarded, as shown in 
Figure 4. The red peak is a multipath of emitter 1, close to the peaks of emitter 2 and stronger than them. 
However, since it is out of the bound defined by the receiver position, this peak can be recognized as a 
multipath and, then, discarded. 

Figure 4. Associating the signals to corresponding receivers. 
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3.2. Range Measurement 

In order to obtain a position fix, several ranges are needed in our system. In the case of additive 
Gaussian noise (AWGN), they can be obtained by the Maximum Likelihood (ML) estimation through 
correlation. Since the pilot subcarriers are modulated by PRBS, the incoming signals can be correlated 
with the local generated replicas, similar to the mechanism used by GNSS receivers [16]. However, 
some differences between GNSS and DVB-T positioning have to be taken into account in the procedure 
of range estimation, due to the different signal structures. 

First of all, the PD has to create a local replica of the SPSs, but the locations of the SPSs of two 
successive OFDM symbols are different, as shown in Figure 2. Therefore, it is necessary to identify the 
SPS location of the current symbol, before creating the local replica. A possible identification method 
is the one proposed in [17] and summarized in Appendix A. This phase can be skipped by adding four 
successive OFDM symbols together [2], as shown in Figure 5. In this case, the correlation can be 
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performed every four OFDM symbols, and the local replica never changes. This can simplify the receiver 
structure and speed up the range measurement in the coarse delay estimate. However, after the coarse 
delay estimate, the correlation should be performed on each OFDM symbol in order to obtain a precise 
delay estimate. 

Figure 5. Scattered pilot location of a quadruplet resulting from combining four 
consecutive symbols. 
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3.2.1. ML Delay Estimate 

It is well known that the ML estimation of the delay of a noisy signal, s r (t) = s(t — T)+w(t) (where 
w(t) is a realization of an AWGN process), is obtained by correlating s r (t) with a signal, s(t), generated 
at the receiver side and by estimating r as the argument of the peak of the correlation function [18]. In 
our case, this rule has to be applied to the signal at the output of the receiver chain shown in Figure 1. 

In the ideal case of a noise-free channel, each received SPS, p, of the k-th OFDM symbol generates, 
after FFT demodulation, a value: 



dt 



c p ae 



(6) 



where a is the attenuation introduced by the channel, and it is a real number; and n is the unknown delay 
to be estimated, normalized with respect to T samp . In the presence of noise, the received value will be: 



d' p = dp + w k 



(7) 



where Wk is the noise contribution. 



To perform the correlation, a local replica of d 1 !, with a variable delay, h: 



P 



c p e 



(8) 



has to be generated and correlated with d' p k . This correlation can be performed in the frequency domain, 
as described in [2,6], as: 



V 



(9) 



P P&Ps(k) 

where N p is the number of SPSs in the OFDM symbol, (•)* denotes a complex conjugate, m = n — n is 
the delay offset between the received signal and the local replica and p s (k) is the index set of the SPSs 
of the £>th OFDM symbol. 
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3.2.2. Correlation in the Ideal Case 

In the ideal case of a single emitter in a noise-free channel, without multipath and frequency offset, 
the correlation, R M (m), can be analytically evaluated for m G A m , where A m is a discrete-time interval 
containing all the possible delays between the incoming signal and the local replica (expressed in number 
of samples). By substituting d'^ = dJl in Equation (9), it is possible to show [9] that the absolute value 

of R M {m) is: 



in / \ I 16 2 

\RM{m)\ = — a a 
y 



sin(7ri?m) 



sin ■ nBm 



(10) 



where is the power of the data and TPS subcarriers and B is defined as: 



iv FFT 

where Pj is the interval between two adjacent SPSs. Notice that R M {m) is a periodic function with a 
period equal to Nfft/Pi- This period restricts the range of delays that can be estimated. From this point 
of view, the quadruplets of Figure 5 are more convenient as Pj decreases from 12 to three. 

In order to mitigate the side lobe of the correlation function, we use a windowing technique. By 
adopting a hamming window, the correlation becomes: 

R W = w E d 'p kp P k Mp) (12) 

r p€p s (k) 

where w(p) = 0.53836 — 0.46164 cos(2np/ (K — 1)) is the hamming window and K is the number of 
non-null subcarriers. 

3.2.3. Delay Estimate in the Case of Multiple Emitters 

In a scenario with more than one emitter, the correlation is characterized by the presence of several 
peaks, each one corresponding to a delay, mi (with / = 1, • • • , N e ), introduced in the correlation by each 
specific emitter. Moreover, each emitter can generate both the line of sight (the only path of interest for 
positioning) and other reflected signals. The peak positions can be evaluated by adopting a threshold 
mechanism applied in the regions of expected delays. The threshold can be calculated as stated in the 
Neyman-Pearson theorem. However, in our simulation, since we are interested in the performance of the 
Position, Velocity and Time (PVT) computation and we know that the signals are present, we skip the 
phase of signal detection, and we simply estimate the propagation delay by considering the position of 
the peak. 

In our experiments, we have used a mechanism consisting of two stages; the first one, called 
acquisition, provides a coarse delay estimate, which will be refined by a tracking stage in the second 
stage. The method adopted in the acquisition stage is based on the comparison of the correlation obtained 
from the measured samples, with the ideal expected correlation, which is, in the ideal case, of no noise, 
and other interfering signals, of the type: 

R Ne (m) = J2 R ( m - m i) (13) 

i=i 
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The criterion adopted in the acquisition stage is the minimization of the mean squared error between 
the measured correlation, RM(m), and the ideal correlation, expressed as: 



where || • || 2 denotes the two-norm of a vector, K is the number of correlation points, m is a test vector, 
m = {m 1 m 2 • • ■ m Nc }, whose elements are variable delays, fhi, introduced for performing the 
minimization procedure indicated in Equation (14), and rh is the vector of the estimated delays. 

The minimization can be done by adopting the Matching Pursuit (MP) algorithm [19], which estimates 
the delays by searching for the correlation peaks in an iterative way. The search space should be the 
whole range, A m , but in our application, the size of the search interval, A m , can be greatly reduced, 
thanks to the initialization phase, which allows a prediction of the peak locations. In this way, the 
estimation of the delays is done only in restricted intervals around the true delays, mi. This allows a 
hot start of the DVB mode, which presents another important advantage: the majority of peaks due to 
multipaths are discarded, as they fall outside the search space identified during the initialization phase. 

Once the coarse delay estimation has been achieved by the acquisition stage, the tracking stage can be 
used to refine this estimation by using an early-late delay lock loop (DLL) similar to the one adopted in 
a GNSS receiver. An example of DLL design can be found in [6,9]. One difference between the tracking 
used here and that in GNSS is the re-acquisition phase. This function can be used to refine the tracking 
performance as the correlation peak jumps from one delay point to the adjacent point. This phase is 
described in Section 3.2.4. 

3.2.4. Re-Acquisition Phase 

Since there is only one sample for each subcarrier, the output of the discriminator will change rapidly 
once the correlation peak jumps from one delay point to the adjacent one. This is different from the 
tracking process used in a GPS receiver, where the receiver can track the signal smoothly from one chip 
to the next one. In order to solve this problem, a re- acquisition function is implemented. If the output of 
the discriminator is above a given threshold, the receiver will enter the re-acquisition stage. This stage 
is similar to the acquisition stage, but the correlation is only calculated upon a small set of delay points, 
which are centered on the delay obtained from the tracking process. Here, we choose three delay points. 
In fact, the receiver trajectory is continuous, and then, the correlation peak can only jump from one 
delay point to one of the adjacent two points (depending on the increase or decrease of the range). If all 
the correlation results are below a given threshold, the tracking is considered as unlocked. The receiver 
will restart from the first step. Otherwise, the receiver goes back to the tracking process with the new 
estimated delay obtained in the re-acquisition stage, as shown in Figure 6. 




2 



(14) 
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Figure 6. The flow diagram of the re-acquisition stage. 
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3.2.5. Final Considerations 

It is well known that the minimum number of reference points for performing the trilateration process 
in a 3-D space is four. In practice, in most applications, three reference points could be enough, as 
the intersection of three spheres gives two solutions, one of which is not realistic and can be discarded. 
Moreover, some aiding data can be provided during the initialization phase, such as the synchronism, 
and the user altitude, further reducing the minimum number of DVB-T emitters. 

In a more realistic situation, three or two DVB-T signals should be sufficient, considering that in most 
applications, some data can be considered quite constant {i.e., altitude) during short/medium GNSS 
outages, especially in a vehicular scenario. Moreover, if the hybrid GNSS/DVB receiver can track its 
position in an area where GNSS satellites randomly appear and disappear, the position provided by 
DVB-T can be used to allow a very hot start of the GNSS-based positioning, even if the GNSS satellites 
are visible in a very short time interval. In general, we can say that a hybrid GNSS/DVB receiver could 
be used for both a hybrid trilateration and for mutual assistance of the two separate positioning systems. 
This paper is devoted to this second aspect, and the simulation experiments with four emitters have been 
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done to prove the feasibility of a positioning method based on synchronous OFDM DVB-T signals, while 
the role of GNSS is only the initialization of the DVB-T-based position computation. 

4. Path Ambiguity Solution 

In this section, two methods are presented to solve the path ambiguity problem, which is caused by 
wrong signal association during the user's motion, which can be observed from Figure 7. One is based 
on pseudorange (PR) comparison and the other on a Doppler aiding decision. These two methods benefit 
from a single visible GNSS satellite. The pseudorange comparison method compares the pseudorange 
estimated from the GNSS signal and the corresponding value computed by using the known satellite 
position and the receiver position in the alternative paths. The Doppler aiding decision exploits the 
Doppler effect in the GNSS signal to assist the receiver velocity and trajectory decision. 

Figure 7. True and estimated trajectories. 
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4.1. Pseudorange Comparison 

When a user receives two or more DVB-T signals with similar or exactly the same propagation delays, 
it can exploit the visibility of a single GNSS satellite to solve the path ambiguity problem by comparing 
the ranges. This method is described by considering the case of two DVB-T signals with the same 
propagation delay. Therefore, two possible positions are calculated. 

In Figure 8, we suppose that the GNSS satellite is at position, S{X s (ti), Y s (ti), Z s (ti)}, with velocity, 
v s , at time tj. The receiver is at position, P 1 {X 1 (t i ), Kx(tj), Zxfa)}, with velocity, v x . Besides this, the 
receiver may also produce a wrong position estimate, which is P 2 {X 2 (ti), Y 2 (U), Z 2 (U)} with velocity, 
v 2 . Therefore, at each time, U, the receiver can calculate the distance between itself and the observed 
satellite. Taking into account the two estimated positions, the two ranges, di and d 2 , are evaluated as: 



di = Vi x s(U) ~ X x (ti)f + (Y s (ti) - Y!(U)) 2 + (Z S (U) - Z^U))* 



(15) 
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and: 

d 2 = V(Xs(ti) - X 2 {ti)Y + (Y 8 (ti) - Y 2 {U)f + (Z S {U) - Z 2 {U)Y (16) 

When the pseudorange comparison method is active, the user can compare d\ and d 2 with p = p + n, 
which is the range derived from the GNSS signal, while n is the error source, which is modeled with a 
Gaussian distribution. It is worth mentioning that the local clock error can be estimated when the user is 
in open sky. Therefore, the range can be obtained instead of the pseudorange. The distance differences 
can be computed as: 

Arfi = d x -p (17) 
Ad 2 = d 2 -p (18) 

Figure 8. Illustration for pseudorange comparison. 
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If one of them is above a given threshold, the corresponding position can be marked as fake, and this 
path will be discarded in the further calculation. We choose three times the standard deviation of the 
satellite range error as a reasonable threshold, considering that the distance difference corresponding to 
the correct path will be less than the threshold, with a probability of 99%. 

4.2. Doppler Aiding Decision 

Another method to select the correct path is based on the Doppler effect, due to the relative motion 
between the satellite and the receiver. The Doppler effect is generally expressed in terms of frequency 
shift, on the basis of the quantities shown in Figure 9. If v r is the true velocity of the receiver, then the 
frequency shift can be written as: 

A/= K-^l|cosfl (19) 
A 
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where \\v s — v r \\ is the magnitude of the relative velocity of the receiver with respect to the satellite, 6 
is the projection angle of the satellite-receiver relative velocity vector to the Line of Sight (LOS) vector 
and A is the nominal wave length of the GNSS signal. 

Figure 9. Doppler effect. 




During the path ambiguity phase, the receiver can compute a vector containing the coordinates of 
position and velocity for each path: one of them is correct, and the others are fake. For each possible 
vector, the receiver can then calculate the corresponding Doppler effect expressed as: 

Af r = H^-^||cos6> r (2Q) 

where v^n-n is the possible velocity for the two positions and 9 re ^ 12 } is the projection angle of the 
satellite-receiver relative velocity vector to the LOS vector for the two possible positions, which can be 
obtained as: 

0 r = arccos ^ — (21) 

\\v s -v r \\- \\P S - P r \\ 

It is known that a GNSS receiver is able to estimate a Doppler shift, A/ GNSS , from the received 
signal. Therefore, by comparing the values of A/ r with A/gnss> the true trajectory can be selected. In 
this method, we do not need to consider the DVB-T receiver clock offset to calculate A/ r , since it is 
directly related to the estimated position. 



5. Simulation Results 



In our application, we assume that the PD is able to receive the signals broadcast by different DVB-T 
emitters belonging to the same SFN. To analyze how the PD can handle these signals, we assume a 
scenario with N e emitters. The distance between the PD and the i-th emitter is denoted as d e ,%, the 
distance, d e ly is associated to the closest emitter and the relative distance, A$ = d e l — d e j, with 
j = 2, ■ • • , N e , is introduced. If the PD is able to receive the closest signal (used as the reference 
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signal), the signals broadcast by emitters with A, < Ap max can be detected, and the relative delays of 
these emitters can be estimated by using the same techniques used for the reception of the useful signal. 
The emitters with Aj > Ap max are seen as noise. In theory, in positioning applications, also, these 
distant emitters could be used, trying to detect them with a dedicated processing. However, they are 
generally blinded by the close emitters; so, in our study, they have not been considered. In our method, 
we assume that the PD starts working in DVB-mode after an initialization phase, which allows the PD 
to evaluate N e and the initial distances, d e ^. 

One important aspect for positioning is the Doppler effect. Since all the emitters are fixed, the 
Doppler effect is introduced by the receiver's motion. In [20], this effect on the DVB-T system has 
been investigated. The results show that for 2 K QPSK, more than 400 km/h is tolerable, but for 8 k 
64 QAM, the tolerance is less than 50 km/h. Therefore, in our application, especially in the 2 K mode 
experiment, the Doppler effect is negligible. 

5.1. Simulation Scenario 

The simulation experiments have been done for the scenario shown in Figure 10, where the user can 
receive four DVB-T signals from four emitters that belong to an SFN, and no GNSS satellites are visible 
when the user enters the grey area indicated in the figure. The user moves along a path with an initial 
location outside the grey area. In our simulation, four DVB-T emitters are placed. 

In the scenario of Figure 10, four DVB-T emitters of the same SFN are placed on the four semi-axes. 
They transmit the same signals on the same frequency, with the parameters summarized in Table 2. 

For the signal generation, we have used the Anritsu MX3700 generator, which is a laboratory Radio 
Frequency (RF) generator that can be modulated by an array of complex baseband samples. As an option, 
a DVB-T/H sample generation software is sold. We used such software to generate DVB-T-compliant 
baseband signals used as simulation sources. The receiver has a constant speed equal to 100 m/s, while 
the direction of the velocity changes from time-to-time. The signal received along the path has been 
simulated by combining shifted and attenuated versions of the generated signal. The delays have been 
simulated by shifting the FFT window according to the corresponding ranges between the user and the 
emitters. The SNR is simulated according to Equation (22): 



where s(t) is the transmitted signal; N 0 is the noise floor; L u is the signal attenuation according to the 
Hata model for urban areas, which is: 



C/N 0 



E[\s(t)\ 2 } 
N 0 L U 



(22) 



L u = 69.55 + 26.16 log(/) - 13.82 \og(h B ) - C H + [44.9 - 6.55 log{h B )} \og{d) 



(23) 



where C H depends on the city size. For small and medium sized cities, it is: 



C H = 0.8 + (1.1 log(/) - 0.7)h M - 1.56 log(/) 



while, for large cities, Ch is: 




8.29(log(1.54/i M )) 2 - 1.1, if 150 <f < 200 
3.2(log(11.75(/i M )) 2 - 4 -97, if 200 < f < 1500 
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where L u is the path loss in urban areas in dB, Kb is the height of emitters in meters, h M is the height 
of the receiver antenna in meters, / is the frequency of transmission in MHz, Ch is the antenna height 
correction factor and d is the distance between the base station and receiver in kilometers. 

Here, we suppose all the DVB-T transmitters are placed on the same height at 100 m. Furthermore, 
since all the emitters use the same frequency, the propagation loss changes only with the distance. The 
SNR of the received signal has been set equal to 10 dB, where the distance between the user and the 
emitter is equal to 5 km. The SNR in the other points is calculated according to Equation (23). 

The position is computed in two dimensions, regardless of the altitude, every 0.2 s. 

Figure 10. Analyzed scenario. 
k Y 




Table 2. Parameters of DVB-T signals. FFT, Fast Fourier Transform; CP, Cyclic Prefix; 
SNR, signal-to-noise ratio. 



Parameters 


Values 


FFT size Nfft 


2,048 


CP length 


1(56 us) 


Signal bandwidth 


8 MHz 


Length of simulation 


20 s, 100 s 


Symbol duration T symb 


280 us 


Sampling period T samp 


109 ns 


SNR 


10 dB 
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5.2. First Experiment 

In the first experiment, the position is evaluated in the absence of multipath and by ignoring the 
propagation loss. The idea is to test only the capability of a system based on DVB-T, and GNSS is only 
used in the initialization phase. The simulation length is set to 20 s. The four DVB-T emitters are placed 
in this way: the two emitters (El (—1,000,0) and E2 (8,000,0)) placed on the x-axis have quite different 
distances with respect to the origin of the reference system, while the other two (E3 (0,4,000) and E4 
(0,-4,000)) on the y-axis are symmetric with respect to the origin. 

Figure 1 1 shows the acquisition results. Four correlation peaks are visible. They correspond to four 
different propagation times with respect to the four DVB- T emitters at the initial location, from where the 
user enters the obstructed area. Since the propagation loss has not been considered, the four correlation 
peaks have the same amplitude. 

Figure 11. Acquisition results (correlation). 
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Figure 7 presents the estimated trajectory. In order to make the figure more visible, the trajectory plots 
are done every 0.4 s, which means that one estimated point between two presented points is dropped. In 
the figure, two different estimated trajectories are visible. The green circles always stay close to the true 
trajectory, while the red diamonds deviate from the true location after the user reaches the origin of the 
reference system. This is because when the receiver reaches the origin, it receives four signals, and two 
of them have the same propagation time (E3 and E4 are symmetric with respect to the origin). Then, two 
different situations may exist in the absence of additional information: 

1. The receiver still makes the correct signal association after this point, and the receiver will 
estimate the trajectory correctly. This is similar to the case of signals from different emitters 
with significantly different propagation times. In this case, the simulation results show good 
performance, even in highly dynamic scenarios (e.g., 100 m/s), as shown by the green circles. 
In our simulation, the receiver achieves a performance with a Root Mean Square (RMS) = 6.8 m, 
Mean error = 6.0747 m and standard deviation a p = 3.0559 m. 

2. The receiver makes a wrong signal association after reaching the origin of the reference frame, 
which means the signal emitted by E3 is associated with E4, and the signal of E4 is considered as 
coming from E3. In this case, the estimated trajectory deviates from the true one, as shown by the 
red diamonds. 

Since we cannot know the true trajectory in the real cases, some additional information (e.g., the cell 
ID information, map information for map matching, signals from other system) is needed to determine 
which trajectory is correct in the case that two or more signals have a similar propagation time. Two 
methods have been introduced in Section 4. 

5.3. Experiment Considering the Hata Model for Urban Areas 

In this experiment, the Hata model introduced in Section 5 is used. Additionally, the path ambiguity 
problem has been solved by the methods introduced in Section 4. In order to simulate a large range of 
SNR, the simulation length is set to 100 s, and the two emitters on the y-axis are located so as to have 
different distances with respect to the origin. 

Figure 12(a) presents the trajectory and the positions of the four emitters, and Figure 12(b) presents 
the correlation results at the initial point of the trajectory in the Hata model. The four correlation peaks 
have different amplitudes, depending on the ranges between the emitters and the PD. From this figure, 
we can see that the amplitude of the correlation peak decreases as the distance between the PD and the 
emitter increases. Therefore, when the distance increases too much, even if it is less than Ap max (the 
maximum detectable range for one DVB-T signal), the ranging operations could become critical if the 
SNR becomes too low. 

By multiplying the estimated delay by the speed of light, we obtain the range: 

r m [n) = r t [n) + e r [n) (24) 

where n is the discrete time normalized by the OFDM symbol duration, r m [n) is the pseudo-range 
measurement, r t [n] is the true pseudo-range and e r [n] is the range error due to noise. The range error 
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for emitter E4 is shown in Figure 13, from which it can be easily observed that the error depends on the 
SNR, which, in turn, depends on the distance between the emitter and the PD. 
In this simulation, the estimated position can be expressed as: 



Px,m[n\ = P x [n} + e x [n] 
PyA n \ = p y[n\ + e y [n\ 



(25) 



where n is the estimation time, Pi[n], with i E {x, y}, is the true value of the PD position on the i-axis, 
Pi,m[n] is the corresponding estimated value and ej[n] is the estimation error. The total positioning error 
can be calculated through the following equation: 



e[n] = yje*[n] + e 2 y [n] 



(26) 



Figure 12. Correlation results in the Hata model. 
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Several different trajectories have been simulated. The true trajectories and the positioning errors, 
e[n], are shown in Figure 14, from which no significant dependence on the SNR can be observed. The 
reason is that along the path, the power of some received signals decreases, while that of others increases, 
so giving approximately a stationary performance in terms of positioning error. 



Figure 13. Pseudo-range error for emitter E4. 




Time (s) 




Sensors 2013, 13 



10213 



Figure 14. Cont. 
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5.4. Experiment with Multipath 

In this section, we simulate the Rayleigh channel for multipath, which is suggested by the DVB-T 
standard [1] for portable reception. Since the method presented in this paper relies on the LOS, we add 
the LOS signal in this model. The related power, delay and phase values can be found in [1]. In order to 
simplify the simulation, we delete the multipath signals, which have a time delay lager than 15 OFDM 
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samples with respect to the LOS signal, corresponding to a 492 m separation in 2 K mode, based on the 
fact that long delay multipaths can be filtered out on the basis of the initial information. 

The correlation with the Rayleigh multipath channel is shown in Figure 15. The Hata propagation 
loss model is also introduced. The trajectory estimation error is shown in Figure 16. From this error plot, 
we see that the position error is comparable with the errors in the absence of multipath. This is probably 
due to the fact that the LOS is present with the strongest signal power. There are four emitters visible, 
which can also improve the performance in difficult environments for 2-D positioning. 

Figure 15. Acquisition results with multipath channel. 
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Figure 16. The position estimation errors with multipath channel. 
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5.5. Experiments with Different SNR 

In this section, we simulate a scenario to test the performance in different SNR. In this scenario, all 
the received signals have the same SNR, which changes from —10 dB to 10 dB for different simulations. 
Additionally, the position errors are shown in Figure 17. In this figure, we can see that when the SNR 
is low, the position errors are slightly worse and the variance is big. When the SNR goes higher, the 
position errors converge to a constant. One thing worth noticing is that the DVB-T signals power is very 
high, above 10 dB, normally due to the error probability requirement of TV service. Additionally, it is 
normal for the DVB-T receiver to receive signals with an SNR higher than 10 dB. 



Figure 17. The position errors in different SNR. 
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5.6. Experiments with One Transmitter Unsynchronized 

In this part, we simulate the scenario in which there is one transmitter unsynchronized. We suppose 
emitter 1 is unsynchronized with the other three emitters, which are synchronized with GPS time during 
the position calculation. Different time offsets are tested; while in each simulation, the time offset stays 
constant. The four trajectories used in Section 5.3 are also simulated here. 

In Figure 18, the mean errors with respect to different time offsets are plotted for the four trajectories. 
We can see that the position error increases, with the time offset increasing for all four trajectories. 

6. Conclusions and Future Work 

In this paper, a positioning method based on DVB-T Single Frequency Networks has been presented. 
This method uses the scattered pilot subcarriers of OFDM symbols to measure the Time of Arrival 
(ToA) of the signals transmitted by the DVB-T emitters. In our study, we suppose that the user is 
equipped with a hybrid GNSS-DVB-T device. GNSS is used in a first initialization phase to solve the 
ambiguities referring to the various DVB-T emitters. When this phase is completed, the user position can 
be obtained by using DVB-T signals when the user enters a GNSS-blocked area. (Note that transmission 
of the emitter ID would be very useful for positioning. Since it is optional, it has not been considered in 
the study, but obviously, it would further simplify initialization and tracking of DVB signals.) 

The method has been tested by simulation in a dynamic scenario. The simulation results show that 
a position mean error = 6.0747 m can be achieved if the user can correctly associate the signals to the 
emitters. If the receiver makes a wrong association, some additional information is needed to determine 
which trajectory is the correct one. The Hata model has been used to simulate emitters with different 
SNRs, and a Rayleigh channel has been also introduced to take into account the multipath. The results 
shows that the range error changes with respect to SNR, while the position estimates do not change 
significantly along the trajectory. However, if one of the DVB-T emitters is unsynchronized, the position 
error will increase as the time offset increases. 

The result of this study is that the structure of the DVB-T SFN signals is such that a PD working with 
DVB-T signals can be conceived. At this point, the performance of the positioning algorithms in the 
presence of different errors sources has to be evaluated. In our future work, a more realistic multipath 
channel will be considered, and an analysis taking into account other error sources will be performed. 
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Appendix 

A. Scattered Pilot Subcarrier Detection 

The detection algorithm proposed in [17] is used to identify the location of the current SPSs. 

This algorithm makes use of the signal property that an OFDM symbol has the same scattered pilot 
subcarriers as the OFDM symbol located four symbols earlier. The correlation of these two OFDM 
symbols is computed, obtaining: 

S(m)= <«~ 4 )*> ™e{3,6,9,12} (27) 

pep s (m) 

where p s (m) is the set of the index of the m-th pattern of the scattered pilot subcarriers and (•)* denotes 
the complex conjugate. 

The maximum value of the four correlations gives the location of the scattered pilot subcarriers of the 
k-th OFDM symbol, which is: 

m — argmax |i>(m)| (28) 

m 

where m 6 {3,6,9,12} denote different scattered subcarrier sets with different first scattered pilot 
subcarrier locations. For example, rh = 3 indicates that the location of the first scattered pilot subcarrier 
is three. The scattered pilot subcarriers of the other OFDM symbols can be obtained by using the rule 
described in [1]. 
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